Superfluid rotation sensor with helical laser trap. 
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The macroscopic quantum states of the dilute bosonic ensemble in helical laser trap at the tem- 
peratures about 10 -6 K are considered in the framework of the Gross-Pitaevskii equation. The 
helical interference pattern is composed of the two counter propagating Laguerre-Gaussian optical 
vortices with opposite orbital angular momenta and this pattern is driven in rotation via angular 
Doppler effect. The macroscopic observables including linear momentum and angular momentum 
. of the atomic cloud are evaluated explicitly. The detection of reference frame rotations using the 

current laboratory equipment is analyzed. 

, PACS numbers: 37.10.Gh 42.50.Tx 67.85.Hj 42.65. Hw 



^ . I. INTRODUCTION 

The superfluids are often considered as an absolute reference frame . For example the circular flows of the liquid 
4 H e in rotating container Q or degenerate quantum gas of a cold bosons in a toroidal trap [H, H| have a quantized 
■ circulation of the velocity along a closed contour L embedded in a superfluid ensemble: 
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where k is a quantum of circulation, <f> is the argument of the macroscopic wavefunction "J/ [5( . When reference frame 
rotates with container with angular velocity the circulation acquires the additional Sagnac term: 

c; 

<j) v ■ dr = <j> (v — £1® x r) ■ df = nn — 2£1® • S, (2) 



where S is area enclosed by contour T. Hence the phase of wavefunction is shifted by Scf>sF = (m/fi.)2f2® • S. Taking 
into account de Broglie wavelength of massive particle moving with speed \v\: Xb = 2ttH/ (m\il\) we get the evident 
similarity of Sagnac shifts for massless photons 5<j)phot and matter waves along with their numerical comparison in 
favour of matter waves: 
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because of smallness of photon's energy hv — hX/c ~ \eV compared to the rest mass of a typical atom. This 
advantage of the matter waves is a somewhat diminished by a substantially larger area ISp/iotl enclosed by optical 
fibers compared to atomic waveguides l^l [6(. Thus studies of atomic interference and search of the novel trapping 
geometries are well motivated. 

The aim of this article is to consider the influence of the reference frame rotations f2® on a superfluid confined by 
helical container aligned along rotation axis z @, H|- The optical trapping by counter propagating Laguerre-Gaussian 
beams (LG) [9( gives the following potential profile for the red-detuned alkali atoms: 
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where the cylindrical coordinates z,r,9 are used, D is the radius of LG [i0| . t is vorticity, k — 2n/X is wavenumber, 
8uj is angular Doppler shift [TTI . [l2j induced by rotation of the reference frame or emulated by rotation of Dove prism 
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in phase conjugated setup (fig.l) [13j. Transformation to the reference frame rotating synchronously with angular 
velocity O = 5ui/2£ = f2$ with trapping helix leads to the time-dependent Gross-Pitaevskii equation (GPE) [5l.l8l.ll4j: 

dty h 2 o 
ih— = - — A* + U opt {z, r, (?) * + g\V\ H - flL z ^, (5) 
Ot 2m 

where the stationary wavefunctions for the superfluid ensemble VP = <fr(z, r, 9) exp{—ijit/K) are given by: 

ti 2 <9$ Airh 2 a<i 

M $ = _ A$ + C^(z,r,d) $ + 2 $ + m^-, g = ^, (6) 
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where as is S*-wave scattering length. We evaluate linear and angular momenta of the superfluid ensemble in a helical 
container Q and discuss the possibilities of rotations detection with this geometry. Hereafter the reference frame 
where laboratory is in rest will be called "laboratory frame" while reference frame collocated with rotating helical 
waveguide where potential U op t is time-independent will be named "observer frame". Noteworthy the "observer" 
velocity vector V with respect to " lab frame" has two components : the azymuthal velocity Vq — fi© x f stands for 
helix rotation around LG propagation axis z while helix pitch velocity V z = {z/ z)8uj/2k is responsible for wavetrain 
translation along z. 



II. TWISTED WAVETRAINS IN ROTATING FRAME 



In order to reveal the basic features of helical confinement let us decompose the potential in the time-dependent 
GPE 0, H| taking into account the paraxiality of laser eigenmodes. Namely the spatial scales in descending order are 
zr » D >> A/2, where A/2 ~ 5 — 0.5/xm, D = 10 — and Zr = D 2 /X (Rayleigh range). Thus the potential 

U op t may be expanded in Taylor series by small parameter z\/(D 2 ) (inverse Fresnel number): 
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This expansion is due to diffractive spread of helical waveguide during motion through LG-beam waist from z — 
point. Consider first the matter wavetrain "J" localized within Rayleigh range z << zr. The atomic waveguide is not 
expanded here, hence GPE is as follows: 
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The variational anzatz r, 9, t) for the order parameter [l6[ in a form of the superposition of the two phase- 

conjugated matter wave vortices [8j] is written for helical trap within LG-beam waist (\z\ << zr) as: 
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where is a longitudinal size of the matter wavetrain, fj, , , are the chemical potentials of the " partial" wavefunctions 
*S?±e. The density p of this atomic wavetrain has a helical shape p(z, r, 9, t) = \^h\ 2 and this helix rotates with angular 
velocity f2 = 2l(fj, f — fJ, b )/h. The qualitative analysis based on anzatz vf^ shows that due to the absence of friction 
the center of mass of superfluid ensemble remains in rest in laboratory frame. But in "observer frame" the center of 
mass will move along straight line z\(t) = zLVbec ■ t = Vit/k ■ t = 5u)/2k ■ t, where Su>/2k is the group velocity of 
the wavetrain (fig.2A). The direction of this rectilinear motion is controlled by a product of the topological charge 
of trapping LG beams I (the handedness of interference pattern) and projection of the vector of angular velocity of 
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"observer frame" fi on z axis. Thus rotating observer will see translational motion of the center of mass of twisted 
wavetrain in the absence of any externally applied force. Noteworthy the order parameter vf^ is a solution of the time 
dependent GPE (eqEJ) which describes purely condensed bosonic ensemble in the absence of thermal background. 

Outside the beam waist the GPE with expanding trapping potential ought to be used (eq[7])-The variational anzatz 
for Y&h > is also a superposition [8]: 
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Because red-detuned solutions are considered, the \^h > is attracted by electrostatic potential U op t towards a beam 
waist \z\ < zr and wavetrain is reflected from low optical intensity regions \z\ > zr towards the center of trap. Thus 
center of mass oscillates with effective equation of motion: 



z 2 (t) +jz-z 2 (*) + [uz/(l + 2ttD /X)} 2 z 2 (t) = F cos(nt) , 
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where oj z = h/(mzR 2 ) is a frequency of oscillations in a "longitudinal" well, F is the averaged amplitude of "kicks", 
caused by reflections from trap boundaries located near z ~ zr, is phenomenological damping constant. The 
dissipation which is absent in initial GPE (eqJT]) appears in this model fcq |ll[) as phenomenological constant, arising 
because of leakage of atoms through potential barrier. The observer rotating with trap will see the onset of wavetrain 
oscillations (fig. 2b) without any externally applied force. Consequently the helical waveguide for neutral atoms will 
transform rotation of lab frame into translational motion of the helical matter wave. 



III. ANGULAR MOMENTUM AND TRANSLATIONS OF SUPERFLUID IN ROTATING HELICAL 

PIPE 

It is instructive to evaluate exactly the macroscopic observables using variational wavefunction > feqs l9U0l) . 
Because GPE (eqsJTHH]) gives coherent mean field wavefunction ^, the expectation values of linear < P z > and angular 
momenta < L z > may be obtained by averaging over macroscopic wavefunction 18]: 

< P z > h =< & h \-ih—\y h >= NH(kf - k b ), (12) 

oz 

taking into account the mutual orientation of the optical angular momenta of LG beams in a trap 

< L z > h =< $ h | - ifi^l*/. >= N£h{l T 1), (13) 
where N is the total number of bosons in a trap, the upper sign stands for the phase-conjugated LG-beams, i.e. for 



helical trap 17|, while the bottom sign stands for sequence of toroidal traps aligned along the propagation axis z. 
As it easily seen from feq fT2l) the ensemble moves along z with momentum h(kf — kb) per particle hence velocity of 
matter wave translation recorded by observer is VL@/{kf — kb). The angular momentum proves to be zero not only in 
observer frame but in a lab frame as well due to disappearance of the moment of inertia for rotating superfluid [5| . 

There exists a remarkable difference between proposed configuration (fig.l) and previously reported Sagnac matter 
wave interferometers @ based on toroidal traps, matter wave grating interferometers (l9j and rotating low-dimensional 
traps where circulation of velocities around closed contour T is essential. The circulation is proportional to the angular 
momentum of a classical particle. In a simplest case of circular rotation in a plane perpendicular to z axis the angular 
momentum is equal to it's projection L z — m\v\r. The quantization of the angular momentum means that operator L z 
averaged over the eigenstates \E' n has discrete spectrum < — ihd/d9\^ n >= ±nh [18]. The quantum ensemble in 
a state ^ with zero circulation or zero angular momentum n — o (a so-called Landau state ) does not feel the rotation 
of container and may be referred to as an absolute reference frame. This behavior had been confirmed experimentally 
both for the liquid helium [l[ and for the trapped alkali gases [20| . 

Nevertheless zero angular momentum state becomes thermodynamically unstable when the speed of angular motion 
fi® of container is increased above a certain critical value [j| . The angular momentum transfer from rotating environ- 
ment to superfliud is facilitated by asymmetrical form of container |l[ or externally imposed ponderomotive optical 
lattice potential [21( but the actual underlying mechanism includes interaction of quantum (superfliud) and classical 
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(normal) components of the ensemble [3|, [5| . This interaction leads to formation of unbounded vortices and vortex 
lattices. For sufficiently fast rotations when fi^ is comparable with classical transverse oscillation frewuency the 
lowest Landau level (LLL) state appears. Noteworthy the vortex lattices also demonstrate a certain independence 
from external rotations. 12211 . 



IV. ROTATION OF LAB ENVIRONMENT DUE TO ANGULAR DOPPLER EFFECT 



Experimentally the rotation of the optical interference pattern Il3l 12(1 has a definite instrumental advantages over 
mechanical rotation The essential features of helical tweezer [23j| are due to apparent topological difference of 
toroidal and helicoidal interference patterns [7](fig.l). The toroidal pattern appears when colliding photons have a 
parallel OAM's while helicoidal pattern [T3| arises for antiparallel OAM's. This follows from direct calculation of OAM 
for the £'th order LG. The OAM is the expectation value of the angular momentum operator L z = —ih [fx V] = 
—ihjfn inside the interaction volume V = irR 2 Z r : 



< L z >(/,(,)=< ¥ (m \L z \*f m >~ J (E+ m )*(-ih[rx V]E+ m ) d 3 f = 

J K. b )r(- ih wo E tf^ rdr - d6dz - eh j^ c > ( 14 ) 

where I(f^b) — e c |E(/,&)| 2 i s the light intensity, b } = \/2e Ey b -, (z, r, 6, t) is the macroscopic wavefunction of a 
single light photon inside a forward (/) or backward (b) LG wavetrains with the winding number t. The retroreflector 
is placed apart from the interaction volume V located near z = plane within Rayleigh range Z T < k^^Do 2 and 
the fields E^ ^ are both calculated inside this volume (fig©. The square modulus ^(/b)! 2 is a probability density 
for energy rather than for a particle (photon) location. This is a Sipe's single-photon wavefunction [24| which is 
proportional to the positive frequency part E^ b of the counter-propagating fields. The usage of the nonrelativistic 

operators of the linear momentum (p = —ihX7) and angular momentum (L z = —ih-J^) is justified by the applicability 
of the paraxial approximation to the slowly diverging LG beams. In this particular paraxial case the spin-orbit 
coupling [25[ may be neglected and the angular momentum of the photon is exactly decoupled to the spin and the 
orbital component: J = S + L. The linear momentum expectation values < p >(/,{,) for the forward and backward 
LG's respectively are obtained in following way: 



<P Z >(/,&)=< *e\P\*e >~ J (E+ m r(-iHVE+ m ) d 3 r = 

It follows from Maxwell equations that the expression for the classical momentum density of electromagnetic field is 
the polar vector P, while angular momentum density M is the axial vector Q: 

P(z, r, 9, t) = e [E x B]; M{z, r, 9, t) = e [ f X [E X B]] (16) 



P, 
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= ±I( fib yirR 2 -Z r /c 2 , 



(17) 



where 2Z r is the length of trapped cloud, R = Dq/2 is the radius if helix. In a same way the classical angular 
momenta L z of LG beams inside the interaction volume V = 2irR 2 Z r projected on z-axis are: 



L. 
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rdr I [ f x [E x B]]dz 



-z r 



± ^' / (/,fc)' 7ri? ■ Z r/(cU(f,b)) 



(18) 



The ratio of the angular and linear momenta proved to be equal to L z /P z <~ I c/uj^f ^ Consequently both classically 
and quantum mechanically when linearly polarized (with (S — 0)) LG-beam is retroreflected from conventional mirror, 
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the polar vector of the optical momentum P changes direction to the opposite one while the axial vector of the angular 
momentum J — L remains unchanged due to the isotropy of reflector 17]. Thus the topological charge I in PCM 
alters sign. 

In contrast to the conventional mirrors the wavefront reversing (PC) mirrors are essentially anisotropic optical 
elements [HI, E3, HE H3 • Due to the internal helical structures, e.g. acoustical vortices in SBS mirrors, the PC mirror 
alters the angular momentum of the incident beam. Thus because the linear momentum P is altered too, the mutual 
orientation of P and L is conserved and the topological charge £ is not changed by an ideal PC mirror. This may 
be also a linear optical loop setup |2 8| provided the even number of reflections per a single loop roundtrip or the 
nonlinear optical PC retroreflector (171 123|| . 

The alternation of the orbital angular momentum in PC mirror is accompanied by conservation identities [l3j for 
the frequencies of the incident and transmitted photons and the projection of the angular momenta L z on z — axis. 
When Dove prism rotates with the angular velocity ft we have: 

I zz • ft + L z = I zz • ft' + L z 

^ / + ^! = ^' + ^, (19) 

where left hand sides of this system correspond to the incident photon and the right hand sides correspond to the 
transmitted one, I zz is the moment of inertia around z — axis. The difference of the angular velocities of the prism 
before and after the photon passage is due to reversal of OAM projection from L z = ±£h to L z = ^f£h ea.([l9| : 

if ■ h 

n-n = _£fL_^. (20) 

From this follows that co-rotation increases the angular velocity of prism. Indeed the energy is transferred to the 
prism by the optical torque \T\ = 2£ ■ P/uif, where P is total power carried by LG As a result Doppler frequency 
shift for the photon u) — LOf is negative in this case: 

8w = J - ujf = ^(ft - ft')(ft' +n) = -2£-si- ^— ^. (21) 

in l zz 

The else consequence of angular momentum conservation is that reversal of angular velocity of the prism ft alters the 
angular Doppler shift sign 5u> and helical trapping potential U „t, (z, r, 9, t) rotates in opposite direction. The remarkable 
feature of this phase-conjugating laser interferometer technique |13l.l30l] is that angular velocity of interference pattern 
ft^ = 8uj/2£ alters when vorticity of LG beam £ is altered (fig. 1) . 



V. CONCLUSIONS 



The properties of superfluid ensemble in helical trap placed in rotating environment were analysed. For ultracold 
bosons (10 _6 K) trapped in vacuum by spiralling optical interference pattern the Gross-Pitaevskii mean field theory 
was applied. The variational approach indicates the possibility of detection of the slow rotations of observer by means 
of measurement of mechanical momentum of atomic cloud. Namely the rotating observer will detect the appearence of 
the translational rectilinear motion and oscillations of matter wavetrain near equilibrium position of the trap without 
any externally applied force. The predicted effect is due to macroscopic quantum interference between two partial 
matter waves with conjugated phase profiles Qj], H3|. The orbital angular momentum of trapping beams [TJ IHJ 
and rotational Doppler shift (l2l . [l3| are the key components for the proposed experiment. 

The GPE-based model considered above is qualitavely applicable to the liquid He 4 flow at (4°K temperature ) 
through a /zm-size capillary [32-34]. In particular the /zm-size helical pipe filled by He 4 superfluid and placed to 
rotating container is expected to eject quantum liquid outwards the pipe with translational speed V = Q,£/k defined 
by angular frequency of container ft, the pitch of the helix A = k/2ir and winding number £. 



[1] G. Hess and W. Fairbank,Phys.Rev.Lett., 19, 216 (1967). 

[2] R.P.Feynman, "Statistical mechanics", Ch. 11, (1972) Reading, Massachusetts. 

[3] A.L. Fetter, Rev.Mod.Phys. 81, 647 (2009). 



6 




FIG. 1: Helical optical trapping potential U op t formed by two phase-conjugated optical vortices with opposite orbital angular 
momenta ±M. Two mutual orientations of the angular momenta of vortices give the interference patterns of right handedness 
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frequency £2 is induced by rotating Dove prism located between trapping volume and phase-conjugating mirror PCM. This 
rotation is equivalent to rotation of the reference frame. The twisted BEC wavetrain in rotating helical atomic waveguide moves 
with a speed Vbec = 0.1/k in a positive or a negative z direction. The direction of motion is controlled by Q.I product. 
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Z axis 

FIG. 2: Twisted matter wavetrain with Gaussian ^-envelope in rotating helical atomic waveguide, a) In a waveguide with 
homegeneous (z-independent) width of potential U op t the wavetrain will move with a constant speed Vbec = ±fM/fc either 
positive or negative direction of z. b) In a waveguide with parabolic z-dependent potential V op t having width Do = £>o(l + 
z 2 /z%}) the wavetrain will experience oscillations between reflection points [29| . 



